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Abstract 

Appropriate combinations of laser beams can be used to trap and manipulate small 
particles with “optical tweezers” as well as to induce significant “optical binding” 
forces between particles. These interaction forces are usually strongly anisotropic 
depending on the interference landscape of the external fields. This is in contrast 
with the familiar isotropic, translationally invariant, van der Waals and, in general, 
Casimir-Lifshitz interactions between neutral bodies arising from random electromag¬ 
netic waves generated by equilibrium quantum and thermal fluctuations. Here we 
show, both theoretically and experimentally, that dispersion forces between small col¬ 
loidal particles can also be induced and controlled using artificially created fluctuating 
light fields. Using optical tweezers as gauge, we present experimental evidence for the 
predicted isotropic attractive interactions between dielectric microspheres induced by 
laser-generated, random light fields. These light induced interactions open a path 
towards the control of translationally invariant interactions with tuneable strength 
and range in colloidal systems. 
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I. INTRODUCTION 


The familiar isotropic dispersion forces between neutral objects arise from random electro¬ 
magnetic waves generated by equilibrium quantum and thermal fluctuations^^^. Depending 
on the context, these forces are known as non-retarded van der Waals-London, Casimir- 
Lifhsitz and, more generally, Casimir forces^^*^. The interplay between Casimir forces and 
electrical double layer forces, which forms the basis of the famous DLVO theory^ describing 
the forces between charged surfaces in a liquid medium, plays a key role in the colloidal 
behavior observed in biological fluids (e.g. proteins, biopolymers, blood cells), foodstuffs 
(e.g. dairy, thickeners, emulsions and creams) or suspensions (e.g. pharmaceuticals, slur¬ 
ries, paints, inks)^. Colloids have also been shown to be extremely well suited for the study of 
phenomena such as crystallisation, the glass transition, fractal aggregation and solid-liquid 
coexistence®^®. External control of isotropic interactions in colloidal systems is therefore of 
key importance. Temperature sensitive swelling of ’smart’ microgel particles offers control 
over soft repulsive forces but the process is slow, shows hysteresis® and the properties of the 
colloids are altered while swelling. In some cases magnetic and dielectric dipolar forces can 
be induced by external helds but these interactions are strongly anisotropic leading to the 
formation of chains or anisotropic domains^®. Other ways to control colloidal interactions 
usually involve the change of composition: Adding and removing electrolytes the range of 
electrostatic repulsions can be tuned and, by dissolving macromolecules of appropriate size, 
attractive depletion forces can be induced^’®. However, despite their widespread and suc¬ 
cessful use, these strategies are still tedious and slow and do not provide the level of control 
over interaction forces that, as discussed here, could be achieved by using external laser 
helds. 

Intense light helds can be used to trap and manipulate small particles^^^^® as well as 
to induce signihcant optical binding forces^^’^^’^® which, in general, are not translationally 
invariant, showing a strong anisotropy that depends on the interference landscape of the 
external helds^®. Here we show that artihcially generated random helds with appropriate 
spectral distribution can provide control over attractive and repulsive isotropic (and transla¬ 
tionally invariant) interactions with tuneable strength and range. In contrast with Casimir 
interactions, where the forces are dominated by the material’s response at low frequencies, 
our results open a new way to explore the peculiar optical dispersion of small particles and 
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artificial metamaterials by selecting the spectral range of the random held. As an example, 
we predict that the interactions between semiconductor particles with relatively high refrac¬ 
tive index, can be tuned from attractive to strongly repulsive when the external frequency 
is tuned near the hrst magnetic Mie-resonance^®. 

Using optical tweezers as gauge, we present experimental evidence for the predicted 
isotropic attractive interactions between dielectric microspheres induced by laser-generated, 
quasi-monochromatic random light helds. We note that isotropic optical forces between 
particles act instantly and can therefore also be applied dynamically. This can potentially 
be useful to anneal defects in periodic structures such as photonic crystals, to increase the 
effective temperature by optically ’shaking’ particles or to stabilize non-equillibrium phases 
such as supercooled liquids and, in general, to control the self-assembly and phase behaviour 
of colloidal particle assemblies on nano- and mesoscopic length scales^’®. 


II. RANDOM LIGHT INDUCED INTERACTION FORCES BETWEEN TWO 
ARBITRARY OBJECTS. 

Early work by Boyer^^ derived Casimir interactions between small polarizable particles 
from classical electrodynamics with a homogeneous and isotropic classical random electro¬ 
magnetic held having the spectral density of quantum blackbody radiation including the 
zero-point radiation held. Here we extend these ideas to “external” artihcial random helds 
with arbitrary spectral density, obtaining an explicit expression for the interactions between 
two arbitrary dielectric objects, which allows a compact description of, random light held 
induced, interaction forces from dipolar (atomic or nanometer-scale) to macroscopic ob¬ 
jects. As a limiting case, when the spectral density of the random held corresponds to 
that of quantum blackbody radiation, we recover the exact trace formulae for Casimir in¬ 
teractions between arbitrary compact objects^®’^®. Related trace expressions have also been 
obtained to describe non-equilibrium Casimir interactions between objects held at diherent 
temperatures^®’^^. 

We hrst analyze the connection between random light induced interaction forces and 
Casimir interactions between two arbitrary objects. We assume that object/particle A is at 
the origin of coordinates and particle B is displaced a distance r along the positive ; 2 -axis 
in an otherwise transparent and non-dispersive homogeneous medium with real refractive 
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index Uh = We consider that the particles are illuminated by a quasi-monochromatic 
random held of frequency u which can be described as a superposition of plane waves 
with random phases and polarizations, propagating in all directions. Each particle can be 
seen as made of discretized, Na and Nb, identical cubic elements of volume v and relative 
permittivity e{u), which act as small polarizable units with an induced dipole proportional 
to the polarizing held, i.e. p(r„,a;) = eoCh Q;(a;)Einc(rn, i^), where a{u) is the polarizability 
given by Q;(a;) = vaQ/[l—ivk^aQ/{Qn)] with 0 : 0 ( 0 ;) = 3(e(a;)—e/i)/(e(a;)+2eft). In the presence 
of a huctuating polarizing held, Einc(r,t), the induced dipoles are huctuating quantities and 
the time averaged force on particle B along the ; 2 -axis may be written as^'^’^^ 



The total force can be seen as the sum of diherent contributions. Although for random 
illumination there is no net force on an isolated particle, the scattered held by B can be 
rehected back by particle A leading to a series of multiple scattering events which give 
rise to a net interaction force between them. An additional contribution arises from the 
correlations between the induced dipoles. At a hrst sight, one could think that the incoming 
exciting helds on the two objects will be completely uncorrelated. However, in a random, 
statistically stationary and homogeneous electromagnetic held, the helds at two distant 
points are correlated, with a cross spectral density of the correlations identical to that of 
blackbody radiation^^. In absence of absorption (when the relative permittivity e(cj) and 
do are real numbers), the sum of the two contributions lead to a conservative interaction 
force which can be expressed in terms of the T-matrix^^ of each individual object [see 
Supplementary Material C] 


F = -VH(r) (2) 

mr) = I I^Tr [in (l - Gb,aTaGa,sTb)] (3) 

where [uE{uj)du] = Ue{oj) is the energy of the huctuating electric held per unit of vol- 
ume and k = nhOjjc is the wave number (c is the speed of light in vacuum). Gb,a is the 
Green tensor connecting the two objects and “Tr” stands for the trace of the T-matrix^"^. 

■k-¥ 

The dependence of the interaction on distance is completely contained in Gb,a whereas all 
the shape and material dependence is contained in the T matrices. The connection with 
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Casimir interactions can be made through Boyer’s approach: When the objects are in equi¬ 
librium with a quantum blackbody radiation, the energy density Ue{oj) corresponds to the 
electric quantum zero-point fluctuations (at zero temperature and positive u) given by^®’^® 
Ue{oj) = UE{oj)duj = hk^ 1 ( 4 : 71 ^)du. For absorbing (emitting) particles, we must include 
an additional contribution to the total force coming from the fluctuating dipoles and the 
corresponding radiated flelds^^ (linked through the fluctuation-dissipation theorem). Inter¬ 
estingly, in equilibrium, this additional contribution conspires with the force due to the held 
fluctuations to give a total interaction potential which is exactly given by Eq. (3), now in¬ 
cluding light absorption and emission (i.e. [Im{e(a;)} > 0) and we recover the exact Casimir 
interaction between arbitrary compact objects^®’^®. In contrast with the traditional Casimir 
forces, equation (3) opens the path towards complete control and tunability of isotropic 
dispersion forces between compact bodies by tailoring the spectral density of artificially 
generated random fields. 

A. Interactions between dipolar electric and magnetic particles 

In the limit of small dipolar particles, Eq. (3) leads to Renne’s result^® obtained from 
quantum-electrodynamic calculations which is identical to Boyer’s^^ based on classical elec¬ 
trodynamics (the more familiar Casimir-Polder result is recovered in the weak scattering 
limit). Recent theoretical works on optical binding (OB) between dipolar particles under 
non-coherent random illumination^®’^° and on dipolar particles near a planar fluctuating 
light source^^ suggested striking similarities between dipolar optical forces in random fields 
and Casimir interactions. However, from a practical point of view, the creation of isotropic 
random light fields and the direct detection of the resulting weak optical binding forces 
between particles at room temperature is challenging. 

For small non-absorbing particles, the interaction energy far from resonance is attractive 
but always much smaller than ksT for realistic power densities, while, at resonance, the 
polarizability is purely imaginary [at resonance, do —t oo, i.e. < 0] which leads to an 
effective repulsion. The latter was shown to play a key role in understanding the collective 
behavior of optical trapped neutral atoms^^. Plasmonic or polaritonic nano-particles show a 
high real and imaginary polarizability close to a resonance but this can lead to a significant 
increase of the temperature^^. Non-absorbing semiconductor nanoparticles, with relatively 
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high refractive index, or colloidal dielectric micron sized particles would offer an attractive 
laboratory to verify our predictions. Semiconductor particles, e.g. silicon spheres with in¬ 
dex of refraction ~ 3.5 and radius ~ 200nm^^, present strong electric and magnetic dipolar 
resonances in telecom and near-infrared frequencies, (i.e. at wavelengths 1.2 — 2/im) with¬ 
out spectral overlap with quadrupolar and higher order resonances^®. Assuming that the 
scattering by these Si particles can be described by just dipolar electric and magnetic fields, 
it is possible to obtain an exact closed expression for the interaction potential, Eq. (3), in 
terms of the electric and magnetic polarizabilities [see Supplementary Material C]. Figure 1 
(I) illustrate how monochromatic random illumination induces a pair potential between two 
identical 230nm radius silicon nanospheres which can be tuned from attractive off-resonance 
(A ~ 2/im ) to strongly repulsive when the external wavelength is tuned near the first Mie’s 
magnetic resonance (A ~ 1.6/im). 

B. Interactions between micron sized dielectric particles 

To derive a simplihed theoretical expression for the interaction potential between micron 
sized dielectric particles we approximate Eq.(3) at lowest order in perturbation theory (see 
Supplementary Material C). In this limit the interaction energy can be seen as given by a 
pairwise interaction wdrs—r^il) summed over the volume of the spheres, i.e. like a Hamaker’s 
integraT, 

U{D) = [ ^u{\rA-rB\) = -y2K{K)U{D,R,K) (4) 

A^a\<R J\rB-re^\<R ^ 

with size independent coefficients K{kn) = 2UE{kn)'nk^ (do/47r)^ as a measure of the mate¬ 
rials contribution to the interaction potential for a given spectral component A;„ = n/jca^/c 
of the random light field. D = r — 2R denotes the gap distance between the surfaces and 
r is the distance between the centers of the two spheres. UE{kn) is the power density 
of the random light field for a specific wave number kn- For monochromatic illumination 
U{D) = —K{k) X U{D,R,k) and in the limit of kR 1 kD we recover the expected 
results for dipolar particles^"^’^® where the interaction energy is proportional to the squared 
volume of the particles. Interestingly, for relatively large particle sizes (in an intermediate 
regime 1 kD kR) we find: 

TT^ R 

U {D) ^ -2UEal — -sm{2kD) , (1 < < kR) (5) 
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FIG. 1: Tunable interactions induced by random light fields, (a) Three representative 
examples of designer random-light-induced pair potentials U{D). (b) Electric field energy densities 
Ue{uj) versus k = UhOj/c. Their corresponding interaction potentials are depicted in (c). (I) Exact 
calculation based on Eq. (3) for nonabsorbing semiconductor particles where the optical response 
is well described by the electric and magnetic polarizabilities. Random light illumination allows 
for switching between strongly attractive (red) and strongly repulsive (blue) interactions. The 
examples show the interaction between two i? = 230 nm silicon spheres (n = 3.5) in water [uh = 
1.33). Based on equation 4 we calculate in (II) the interaction potential due to a monochromatic 
random light illumination between uncharged dielectric particles at different wavelengths (compare 
the corresponding colors). By changing the wavelength (or wave-number k) the range of the 
interaction potential can be tuned precisely. (Ill) Properly designing the energy density spectrum 
of random light helds allows for exponentially attractive dispersion potentials (black curve). This 
property can be used to produce Morse type potentials on a colloidal level using dielectric-charged 
particles. The exponential electrostatic double-layer repulsion calculated for a contact potential of 
180 k^T and a Debye length of 20 nm (blue) and the exponentially attractive dispersion interaction 
(black) superimpose to a colloidal Morse-potential (red). 

i.e. for relatively large particle sizes, our approach predicts that the interaction energy 
should scale with the particle’s radius. Otherwise Equation (4) can be easily computed 
numerically. Figure 1 (11, 111) illustrates the predicted theoretical sensitivity of the induced 
pair interactions between low-index (n=1.68) micron-sized spheres to changes in the wave¬ 
length of the illuminating random held. Since forces can be strongly frequency dependent. 






















actual dispersion forces could be manipulated by tuning the spectrum of the random light 
held. As an example, in Figure 1 (II) we show that it is possible to tailor the range of an 
effective attractive dispersion interaction potential with constant depth by selecting different 
wavelengths and held power densities [calculations are performed using Eq. (4)]. Selecting 
an appropriate continuous spectrum of the random light helds it is also possible to design 
a nearly exponential attractive potential. Taking into account the well known exponential 
electrostatic double layer repulsion, it would theoretically be possible to induce an ehective 
Morse-type interaction potential on a colloidal level, as shown in Fig. 1 (III). 


III. EXPERIMENTAL OBSERVATION OF RANDOM LIGHT INDUCED INTER¬ 
ACTIONS BETWEEN DIELECTRIC MICROSPHERES 

Experimentally, the manipulation of relatively large, micron-sized, dielectric particles 
requires high intensity laser helds: In conventional optical tweezing experiments, the focussed 
laser intensity Iq = P/A required to tightly trap a micron-sized polystyrene sphere (refractive 
index n ~ 1.6), is on the order of mW/pm^ and the corresponding power density is Ue = 
rih/c ■ Iq. Therefore, to create random light helds with a comparable power density the 
surface area A cannot exceed some tens of micrometers squared, for an incident laser power 
P on the scale of Watts. To this end we have designed a miniaturised sample cell that allows 
for the simultaneous creation of a random light held in a small cavity by illumination with 
a strong green laser A = 532nm and the observation of optical binding forces between two 
isolated microspheres (Figure 2). We use the technique of time shared optical tweezers^®^^^ in 
combination with umbrella sampling^® to probe the particle pair interaction potential U (r) 
of two melamine microspheres (n=1.68) as a function of their center-to-center separation 
distance r. The principle of this method is to trap, at the same time, two particles in two 
identical optical tweezers that are separated by a distance tq using a near infrared laser beam 
with a wavelength of 785 nm. By monitoring the relative motion of the trapped particles 
with a digital camera information about the pair potential P(r) around tq is gained. In 
contrast to vacuum, the interactions between particles suspended in water commonly involve 
both van der Waals and screened electrostatic repulsive interactions (DLVO interactions). 
The latter ensures the stability against particle coagulation. Equally, in our measurements, 
this repulsive part dominates at very short distances and in turn this allows us to probe 
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FIG. 2: Experimental setup. An intense green laser (A = 532nm) is weakly focussed on one 
side of the sample cell using a / = 75mm lens. A turbid scattering layer of thickness 20/im at 
the entrance of the cell creates a random field distribution inside a light filled cavity of dimensions 
~ (50/rm)^. From the opposite side the sample cell is illuminated with a tightly focussed near 
infrared laser (A = 785nm) creating a set of two time shared optical traps. Both laser beams can 
be steered by galvano mirrors. A white light source is employed for the broadband illumination 
of the cell allowing the tracking of the particles by video microscopy. Different filters are used to 
spectrally isolate the different optical paths. The interaction potential between the two trapped 
spheres is obtained by monitoring the thermal Brownian motion inside the optical traps with a 
digital camera (CCD). The inset shows an enlarged view of the sample cell: Two micron-sized 
colloidal particles suspended in water are trapped in the center of water filled layer of thickness 
t = 15^m. The water layer and the turbid layer are separated by a glass wall of thickness h ~ 20/im. 
The dichroic mirror at the bottom of the clear layer reflects more than 99% of the incoming light 
leading to multiple reflections inside the cavity. 

the superimposed light induced attractions, without particles sticking together irreversibly. 
We follow the approach of Grier and coworkers^® to obtain an autocalibrated measurement 
of the colloidal interaction potential by analysing the differences between the distributions 
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of the particle positions in the presence and absence of optically induced forces (see also 
Supplementary Material A and B). Thus in our experiments, by turning on and off the 
random light field, we obtain directly the light-induced interaction potential without an 
explicit measurement of the complete pair interaction potential. Because of the hnite width 
of the laser traps the particles sample only a small range of U around tq. In order to 
determine the potential over a wider range of distances we repeat the measurements for 
different relative positions of the optical traps (umbrella sampling). Typically we start at 
an average trap separation that is near contact and increase tq in steps of A = 40 nm until 
6 different trap positions are scanned. A is chosen to be smaller than the width of the trap 
potential 2a ~ 140 nm to have sufficient overlap between neighbouring trap separations. 

The final results for U{D = r — 2R), shown in Figure 3, are obtained by averaging over 
16 independent experimental runs carried out under the same conditions. We hnd clear 
evidence for optical binding in a random light held. Moreover, we are able to quantify the 
potential U{D = r — 2R) as well as its dependence on the incident light power. Varying 
the laser power from 5.0 W to 1.0 W the attractive potential weakens. For an estimate 
of the contact potential Uo = U{D = 0) we adjust the prefactor in Eq. 4 for a best ht 
to the data as shown in Figure 3. The overall agreement between experiments and theory 
is remarkable except for the predicted oscillation of the interaction potential at large inter 
particle distances. This could be partially associated to the absence of a full dvr isotropic 
illumination in the experiment due to the absence of incoming photons with momentum 
near parallel to the mirror. It could also be due to the statistical average over different 
experimental realizations as positions extracted from these different measurements can be 
slightly shifted with respect to each other, which may result in a smearing of the curves. 
Finally we note that the slight mismatch with theory could also be due to the approximations 
made when deriving Eq. 4. 

In Fig. 4 we display results for the contact potential Uq as a function of laser power P 
for the three different particle sizes. From the slope of the linear hts (dashed lines) we can 
extract the normalized contact potential Uq/P which does not depend on the laser power. 
Uq/P increases with particle size and the data set is consistent with the linear increase as 
predicted by Eqns 4 and 5. Finally we attempt a quantitative comparison between the 
experimental results and the theoretical predictions for the contact potential Uq. We note 
that such a comparison is based on a number of uncertainties related to the approximations 
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FIG. 3: Optical binding between dielectric microspheres in random light fields. From the 
analysis of the the thermal motion in the two adjacent optical traps the particle interaction potential 
U{D) is obtained experimentally. The melamine (n=1.68) microspheres have sizes 2R = 2 /rm (a), 
2R = 3 /im (b) and 2R = 4 fim (c). All particles are suspended in a 2.7 mM KCl aqueous solution. 
In the figure the measured interaction potential is shown for three different laser power settings 
P in each panel. The vertical lines show the estimated contact positions and the experimental 
uncertainty AD ~ ±5nm. By numerical evaluation of Eq. 4 we can fit the data and then extrapolate 
U{D;P) to obtain values for the contact potential Uq as a function of particle size 2R and laser 
power P. The inset shows an enlarged plot for 2R = 4/rm and P = 5W. 

made in the theory as well as the experiment. From Eq. 4 we compute the theoretically 
predicted values for melamine microspheres with an index of refraction of np = 1.68 in 
water {jih = 1.33) as a function of particle size. For the cross-sectional area of the light hlled 
cavity we use A ~ 0.004 mm^. Moreover we take the light power in the cavity equal to the 
incident laser power for a vacuum-wavelength Aq = 532 nm. Under these assumptions the 
theory predicts jp ^ o.2 ■ kpT/W x (2i?//im). This results matches the experimental 
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FIG. 4: Laser power and particle size dependence of the contact potential. Contact 
potential Uq as a function of the laser power P for particles of different size 2R = 2/nn, 3;um, 4/im. 
From a linear fit (dashed lines) we estimate the normalized contact potential Uq/P for the respective 
particle sizes. The inset shows the dependence oiUo/P on the particle diameter 2R and the dotted 
line is a linear fit to the data t/o/P ~ 0.17 • ksT/W x {2R/jim). 


value Ip 0.17 ■ ksT/W x {2R/iim), inset Figure 4. Given the approximations 

made such a near quantitative agreement might be somewhat fortuitous but nonetheless the 
overall agreement between experiment and theory supports our hndings. 
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Appendix A: Experimental Methods 


A focused green laser beam (Coherent Verdi-V5, A = 532 nm, calculated 1/e^ beam- 
waist 2wo = 25 /im) with a laser power of up to 5 Watts is focused with a / = 75 mm 
lens on the back surface of a glass capillary hlled with a dense amorphous solid composed 
of PMMA (Polymethylmethacrylate) beads, diameter ~ 0.4/im, with a layer thickness of 
~ 20 fim. This hrst layer scatters more than 99 % of the incident power and thereby creates 
a random light held in the sample cavity of thickness t ~ 15 pm (inset, Figure 2). The latter 
is hlled with a very dilute suspension of melamine microspheres (Microparticles GmbH, 
Germany) with a diameter of 2 pm < 2R < 4 pm and refractive index n = 1.68, dispersed 
in aqueous solution and sealed with UV-curable glue. In order to screen the electrostatic 
repulsions we add 2.7 mM KCl. This results in a Debye screening length of Ad — 6 nm. 
The thickness of the layer containing the melamine microspheres is controlled by adding 
a small amount of t = 2R = 15 pm silica spheres that serve as spacers. Both layers are 
separated by the outer glass wall of the capillary with a thickness of h = 20 ± 2 pm (GM 
Scientihc). A dichroic mirror with a rehectivity > 99% at A = 532nm (for incident angles 
0° to 45°, custom made by Asphericon GmbH, Jena, Germany) is placed on the opposite 
side of the cell. The mirror is transparent for wavelengths above A ~ 650 nm allowing 
both the trapping and the visual observation of the melamine microspheres. The clear 
layer sandwiched between the turbid layer and the mirror (cavity) thus has a thickness of 
h -|-1 ~ 35 pm. The thin slab geometry together with the scattering in the turbid layer and 
multiple reflections in the cavity assure that in the center of the sample cell we generate 
a fairly isotropic and quasi-monochromatic random light held. The volume speckle inside 
the cavity displays intensity huctuations on a typical length scale^^ of A/2 ~ 200nm. In 
contrast with the anisotropic static light speckle pattern used to generate random potentials 
landscapes in the quest for Anderson localization of matter waves"^^, we generate random 
temporal huctuations of the light helds by rapidly scanning the green laser focal spot over 
the surface of the turbid layer. We set the scan distance to p ± 7 pm at at an oscillation 
frequency of 500 Hz using a galvano mirror. The oscillation frequency is chosen in a way that 
light held huctuates randomly much faster than tb = R^/{6Dq) = 776 ms, the time scale 
for Brownian motion of the smallest melamine microspheres under study. Here Dq denotes 
the Brownian dihusion coefficient. It has been recently reported that the collective motion 
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of a large set of microspheres under the influence of both Brownian and self-interacting 
optical forces becomes active and their dynamical quantities are no longer representative 
of thermodynamic eqnilibrium^^. In our case, however, the rapidly fluctuating speckle is 
not coupled to the motion of the melamine spheres. It is important to note that roughly 
2/3 of the incident laser light will be reflected by the turbid layer. However, the reduced 
incident power is compensated by mnltiple reflections inside the cavity except for the small 
residual losses by the dichroic mirror. Therefore we expect the laser power in the cavity to 
be approximately the same as the incident laser power P. Reflections in the cavity can lead 
to an increase of the effective surface area A as the reflected light can spread ont laterally. 
To obtain an estimate of A we image the residnal green light, for an incident laser power of 1 
W, transmitted by the dichroic mirror by increasing the exposure time and the gain setting 
of the camera corresponding to an increase in detection efficiency by more than fonr orders 
of magnitude. From an analysis of the slightly elliptical intensity distribution we derive an 
areal cross section of H ~ tt x 41 x 34yum^ ~ 0.004mm^ based on the 1/e decay length along 
the major and the minor axis. This means that for a laser power of up to 5 W we can indeed 
reach intensities ~mW/comparable to the case of optical tweezing. 

In one experimental run we observe the Brownian motion of two micron sized melamine 
microspheres that are held at a mean distance tq in the center of the light hlled cavity, both 
laterally and axially, using a specihcally adapted Nikon Eclipse TSIOO bright held microscope 
composed of (i) a long working distance objective (20x/0.42 EO Plan Apo ELWD) for the 
sidewise white light illumination (ii) an oil immersion objective for both the trapping and 
the observation of the particle motion using a CCD camera and (iii) a notch hlter to hlter 
out residual stray light coming from the trapping laser as well as a dielectric mirror to conple 
the trapping laser into the optical path of the microscope. Note, the bright held illnmination 
takes place across the hrst dihnsing layer. Eor the laser trapping we use a near infrared laser 
beam (Toptica DL 100) with a wavelength of 785 nm and a power of ~8 mW (measured at 
the exit of the laser) that is rapidly switched between two positions with a galvano mirror 
to produce the time shared dual-traps^^^^^. 

Using a digital camera (Prosilica GC650, Allied Vision Technologies GmbH, Germany) 
we record images of 120 x 120 pixels with a frame rate of 90 Hz and an exposure time of 
0.3 ms at 80X magnihcation. The edge length per pixel is dpix ~ 100 nm which provides a 
snb-pixel localization accnracy of the particle center better than lOnm^®. For a given mean 
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separation distance ro we perform two experiments (see Figure 2). In a first experiment we 
acquire a movie of 4000 images at a frame rate of 90 Hz under the influence of a random 
light held. In a subsequent reference experiment the random light held is turned oh and the 
measurement is repeated under otherwise identical conditions. The recorded sequence of 
images is analysed using a standard particle tracking algorithm^^ to obtain distributions of 
the center to center separations of the trapped particles for both experiments. We follow the 
approach of Grier and coworkers^® to obtain an autocalibrated measurement of the colloidal 
interaction potential by analysing the diherences between the distributions in the presence 
and absence of optically induced forces. 

Appendix B: Supplementary Material: Experiment 

1. Laser trapping experiment and data treatment 

The focused beam of a Topica DL 100 diode laser operating at a wavelength of 785 nm 
is time shared between two points using a galvano mirror (Galvoline G1432) driven by a 
square-wave oscillation at a frequency of 500 Hz. A telescope (Thorlabs, 3x Galilean optical 
beam expander, BE03M-B) is used to match the beam size with the back aperture of an oil 
immersion objective (Nikon 60x PlanApoVG, N.A.= 1.4)^^. The location of the telescope 
is chosen in a way that the back focal plane of the oil immersion objective is imaged onto 
the galvano mirror, which allows for identical dual-traps^®. Finally the time shared beam 
is focussed into the water layer of the sample cell to form the dual-trap. The particles are 
trapped in the middle of the water layer to minimize wall effects. The average distance 
between traps’ centers can be changed by adjusting the amplitude of the galvano mirror 
oscillations. With the GGD camera we record images of 120 x 120 pixels with a frame 
rate of 90 Hz and an exposure time of 0.3 ms. With a micro-scale the effective pixel size 
is measured to be dpix ~ 0.1 /rm. The recorded images are analysed using an adapted 
MATLAB (The MathWorks, Inc., USA) code based on the particle tracking algorithm by 
Grocker and Grier"^® to hnally obtain the particle positions on each picture. We quantify the 
transversal instrumental resolution of our apparatus by tracking two 2R = 2 /xm particles 
with a center-to-center separation of that are permanently adsorbed to the lower glass 
surface of the water layer. A Gaussian fit to the measured distribution reveals a standard 
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FIG. 5: Thermal motion inside adjacent optical traps. Two particles with radius R inside 
two identical optical traps are positioned at a distance tq. For an isolated particle the distribution 
of positions is Gaussian with a standard deviation a, which is set by the trap stiffness. Thus 2a is a 
measure for the typical distances r — tq probed by thermal motion. Interactions between particles 
lead to a characteristic change in the distribution of particle positions. Attractive interactions 
increase the probability for the particles to approach. Precise measurements of fpairif) are therefore 
a sensitive tool to determine the particle-particle interaction potential U{D = r — 2R). 


deviation of axy = 0.077 pixel = 7.7 nm reflecting the transversal instrumental resolution. 
Moreover we have verified that out-of-plane fluctuations due to the finite trapping strength 
are negligible in our experiment'^®. 

The thermal motion of the two particles in the adjacent traps is illustrated in Figure 
5. For a given mean separation distance tq of the time shared optical traps we perform 
two experiments (see Figure 1 of main text). In a first experiment we acquire a movie of 
4000 images at a frame rate of 90 Hz under the influence of a random light field. In a 
subsequent reference experiment the random light held is turned off and the measurement 
is repeated under otherwise identical conditions. The recorded sequence of images is anal¬ 
ysed using a standard particle tracking algorithm^® to obtain distributions of the center to 
center separations of the trapped particles for both experiments. We follow the approach of 
Grier and coworkers^® to obtain an autocalibrated measurement of the colloidal interaction 
potential by analysing the differences between the distributions in the presence and absence 
of optically induced forces. 
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where f^i and /o are the corresponding distribntions of center-to-center separations. 

From n measurements of the center-to-center separation r we compute the pair distribution 
fpairi'f') using the technique of nonparametric density estimation^®’^®: 
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where r* reflects the separation distance determined from one image i (one measurement); 
hopt is a smoothing parameter. The estimator’s kernel J[(r — ri)/hopt] can be any smooth 
function that satishes the following conditions: (i) continuous and symmetric around zero 
(ii) integrable with its maximum Jmax at zero and (iii) normalized and non-negative*^®. For 
convenience we choose a Gaussian function of the form: 
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The smoothing parameter hopt reflects the kernel’s bandwidth. A proper choice of hopt 
is crucial. A too large width obscures features in the pair distribution fpairij') whereas 
a too small width yield noisy results. A good trade-off is given by Silverman’s rule^®: 
hopt = where <7^ is the standard deviation of all separation distances r*. The 

beneht of nonparametric density estimation over histograms is (i) the convergence speed; for 
n data points the statistical error in histograms decreases as whereas for nonparametric 
density estimation the error improves as (see Refs. 39,47). More importantly the 

nonparametric density estimation does not rely on the choice of discrete bins. 


2. Total interaction potential of the charge stabilized microspheres 

Particles suspended in water involve both van der Waals and double layer (dl) electrostatic 
repulsive interactions which in combination can be described by the well known DLVO 
theory^. The latter is dominantly repulsive for stable suspensions and thus prevents particle 
coagulation. Equally, in our measurements, this repulsive part dominates at very short 
distances and in turn this allows us to probe the superimposed light induced attractions, 
without particles sticking together irreversibly. For illustration we show in Figure 6 a typical 
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FIG. 6: Calculated interaction potential for charge stabilized colloidal particles of size 2R = 2/rm 
suspended in water in the presence and the absence of the random light (rl) induced interactions. 
The combination of van der Waals attractions (vdW - red line), Hamaker constant A = 0.1 k^T, 
and electrostatic double layer repulsions (dl-blue line), Debye length Ad = 6 nm, contact potential 
Udi{D = 0) = 44 k^T leads to the DLVO potential (black line). The attractive potential due 
to random light forces (green line), monochromatic random illumination with a contact potential 
of -2 k^T, Eq.(6) is superimposed leading to the total potential (magenta line). The repulsive 
part dominates at very short distances and in turn this allows us to probe the superimposed light 
induced attractions, without particles sticking together irreversibly. 

DLVO potential for micron sized particles {2R = 2/im) and the superimposed attraction due 
to random light helds corresponding to the case shown in Figure 3 (a) in the main text. Exact 
values for the Hamaker constant and the contact potential Udi{D = 0) are not known and 
we have chosen reasonable estimates consistent with the observed stability of the melamine 
particles. 

3. Amorphous turbid layer 

The turbid layer at the entry of the light hlled cavity is composed of a dense amorphous 
assembly of PMMA (Polymethylmethacrylate) particles, diameter ~ 0.4/im. We prepare 
the sample by hlling a hollow rectangle borosilicate glass capillary (CM Scientihc) with a 
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height of ~ 20 /iin and a width of w = 200 /xm with a concentrated colloidal suspension 
with a particle volume fraction of approximately 0 ~ 0.35 and then let the suspension dry. 
We perform scanning electron microscopy on the dry particle layer by breaking the capillary 
after the experiment. The images (not shown) reveal a densely packed random structure in 
the bulk of the dried sample and a thin boundary layer with a crystalline structure close 
to the cell wall. We measure the line-of-sight transmission by collimating the 532nm laser 
beam and masking it with a 50 fim pinhole that we place as close as possible (ca. 1 mm) to 
the front surface of a glass capillary. On the opposite side we record the far held intensity 
prohle by directly placing the sensor of the digital camera. We estimate the line of sight 
transmission to Tios = 0.1 %. We estimate the total transmission by placing a high numerical 
aperture objective (Nikon 60x PlanApoVC, N.A.= 1.4) on the opposite side of the sample 
cell in order to maximise the acceptance angle for transmitted light. The collected light is 
then projected on a screen, imaged and analyzed with the digital camera. From this we 
obtain an estimate for the diffuse total transmission of Td*// ~ 1/3. 


Appendix C: Supplementary Material: Theory 

1. Field-Field correlations and cross-spectral density in a stationary random field 

We consider a huctuating electric held, E(r,t) in a transparent and non-dispersive ho- 


tiotemporal huctuations. 
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where Gij{r, r', u) are the matrix elements of the Green tensor. 
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(I 3 is the identity tensor) and UE{r,u) is related to the time-averaged electric energy per 
unit volume. 


{U{r,t)) = ^(|E(r,t)|2) = UE{uj)du. 


(C4) 
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2. Multiple scattering between two compact bodies 


We consider a particle, A centered at the origin of coordinates and a particle B displaced 
a distance r along the positive ; 2 -axis. From a physical point of view, instead of a continuous 
approach, each particle can be seen as made of discretized, Na and Nb, identical cubic 
elements of volume v. This is also known as a discrete dipole approach (DDA)^®. In the 
presence of an external polarizing held, Einc(r,n;), each volume element acts as an induced 
dipole proportional to the polarizing held, i.e. 


p(r„,a;) = eoSt a(a;)Einc(rn, ca) 


(C5) 


where a{u) is the polarizability which, for cubic or spherical elements of volume v, is given 
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The polarizing held on a given element in particle B, Einc(r^, w), is given by the solution of 
the multiple scattering problem; 


Nb Na 

. n n , . n , ^ 


Einc(rf) = Eo(rf) + ^ G(rf, r®)Ei„c(r^) + G(rf, ri)'Einc{r^) (C7) 

m^n m 

(and an equivalent equation for particle A). For simplicity in the notation, we do not include 
here the explicit ca-dependence. These are a set of 3Na + 3Nb equations that can be written 
in compact matrix form as 


Einc(-B) — Eo(-B) + aBk‘^GB,B^mc{B) + aAk‘^GB,AEimciA) 
Einc(A) = Eo(A) + a Ak"^ G A,A^inc{,A) + aBk'^G A,B^inc{,B) 
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Introducing the T-matrix, dehned as 

T-‘(r® 0 = ^13 -G(r® 0(1-i„„), or s - Gb.b (CIO) 

where Isat^ is the 3Nb x 3Nb identity matrix (and an equivalent expression for T^), the 
formal solution of the scattering problem can be written as 
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Our approach can be seen as the DDA-like version of the well known T-matrix approach 
of multiple scattering of electromagnetic waves by two different objects usually described in 
terms of a basis of multipolar vector wave functions (see for example Ref.^^’^°). 


3. Optical interactions between two compact bodies induced by random light fields 


In the presence of a random (stationary) field, Eo(r, t), both the dipoles and the polarizing 
helds are, in general, fluctuating quantities and the time averaged force along the ; 2 -axis may 
be written as the sum of two different terms [see for example, Ref.^’^] 
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where the first term describes the force induced by the fluctuating (external) field, with 
the corresponding induced dipole as discussed in the main text. The second involves 
the (spontaneous and thermal) fluctuations of the dipole 

We focus on lossless particles and discard the second contribution (in absence of absorp¬ 
tion, there are no spontaneous and thermal fluctuations of the dipoles). From Eqs. (C5) 
and (C13), the total time-averaged force on particle B is then given by 
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where Einc(rn;i^) is given in (C7) and the gradient of the incoming field is the sum of three 
different terms 
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These three terms give three different contributions to the total force on particle B. The 
first term, can be seen as coming from the homogeneous radiation field on particle B 
(which after arriving at R, suffers multiple scattering events with particle A). The second, 
comes from the radiation Erst scattered by particle A. The last term, arising from the 
multiple scattering interactions inside the particle, does not contribute to the total force on 
B since these interactions cancel out when summing over all the dipoles in B after averaging 
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over the random field. Taking into account that 
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we find Ff = du[F^^\u) + F^\u)] with 
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where “Tr” stands for the trace of the 3Nb x 3Nb matrix. After some algebra and taking 
into account that in absence of absorption (i.e. e{u) and do are real) 


ImT^i = ^IsTVs - ImGR,R , ImT \r^,r^) = -ImG(rf,r®), 
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the second contribution can be shown to be given by 
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Adding (C18) and (C20) we finally obtain that, in absence of absorption, the total force is 
conservative F^ = —dU{r)/dr with an interaction potential given by 
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The dependence of the interaction on distance is completely contained in Gb,a whereas all 
the shape and material dependence is contained in the T matrices. 

In the case of equilibrium thermal blackbody radiation the electric energy density, Ue{(^), 
is given by^^’^® 
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which, at zero temperature gives ue{oo) = h/c^/(47r^). For absorbing (emitting) particles in 
equilibrium, we can include in Eq. (C13) the contribution of the fluctuating dipoles and 
the corresponding radiated fields^^ (linked through the fluctuation-dissipation theorem). 
Interestingly, in equilibrium, this additional contribution conspires with the force due to the 
field fluctuations to give a total interaction potential which is exactly given by Eq. (C21), 
now including light absorption and emission (i.e. [Im{e(ci;)} > 0) and we recover the exact 
Casimir interaction between arbitrary compact objects^®’^®. 
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4. Attractive and repulsive interactions between dipolar electric and magnetic 


particles 


Submicron dielectric spheres made of moderate permittivity materials present dipolar 
magnetic and electric responses^®, characterized by their respective hrst-order “Mie” coeffi¬ 
cients, in the near infrared, in such a way that either of them can be selected by choosing 
the illumination wavelength. The scattering properties of Silicon and other semiconductor 
nanoparticles^®, can be well described by their electric and magnetic polarizabilities, being 
negligible the contribution of higher order modes (contribution of higher order modes can 
be relevant when the interparticle distance D becomes of the order of the particle radius -a 
detailed analysis of these interactions will be described elsewhere-). 

When the optical response of the particles can be described by their electric and magnetic 
polarizabilities, a^{uj) and respectively (n = A,B). The presence of an external 

polarizing held induce both electric, p and magnetic, m, dipoles, i.e. 
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where Z = \//io7(d)^ is the impedance of the homogeneous medium. The polarizabilities 
are simply related to the hrst electric, Oi, and magnetic, bi Mie coefficients: 


Gtt 








(C25) 


When we just consider dipolar particles we can write 
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with 
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and g{r) = e*^'’/(47rr) the scalar Green function. 

In absence of absorption the trace formula for the interaction potential [Eq. (C21)] can 
be calculated in closed form as: 
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As shown in Figure 2 in the main text, for two identical particles near the hrst Mie dipolar 
magnetic resonance the interaction force can be repulsive in analogy with the repulsive 
interactions between resonant atoms^^. 


5. Weak scattering approximation 


In the weak scattering limit, we can expand (C21) 
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(C32) 

(C33) 

(C34) 


It is easy to see that this is the result obtained at lowest order in the so-called Born (or 
Rayleigh-Gans-Debye) approximation. For two identical spheres of radius i?, their centers 
being a distance r apart, in a quasi monochromatic random held, the interaction energy in 
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the weak scattering limit can be seen as a Hamaker’s integrah^ 


U{D) = -K xl({D,R,k) 
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(C36) 


with D = r — 2R and K = {du 2uE{(x)}nk^ (Q;o/(47r))^. 

6. Random light forces between dipolar electric particles: gravitational like inter¬ 
actions 


If we consider the long wavelength limit, where the magnetic polarizability is negligible, 
Eq. (C21) takes the simple form 

U{r) = ‘^Ue{u) Im I In ([l - {a^k^Gu{r)Y] ) I (C37) 

Ki=x,y,z J 

which can be shown to be equivalent to Eq. (11) in Ref.^° in absence of absorption. 

A remarkable prediction concerning optically induced interactions between atoms, 
molecules or small dipolar particles^^ is that, after averaging over all orientations of the 
inter-atomic axis with respect to the incident beam, the interaction is an isotropic long- 
range, “gravitational-like”, 1/r potential in the near held. It was suggested^^ that this 
averaging could be experimentally achieved by an isotropic external illumination by means 
of multiple incoherent beams which, for atomic systems, could give rise to stable Bose- 
Einstein condensates with unique static properties'^. An alternative is to average over all 
orientations and polarizations of the incoming, uncorrelated, plane waves^^; In the weak 
scattering limit, expanding Eq. (C37) leads to 


27r 


77(r) RiIm <j ^ (a''k^Gn(r))‘ 


(C38) 


, i=x,y,z 


which in the short distance limit gives the above mentioned long-range 1/r dependence of 
the optical interaction potential in agreement with previous results^^’^^. It is worth noticing 
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that similar ideas were considered in the earlier proposal by Spitzer^^ of the so-called mock 
gravity, gravity-like interactions between matter in the universe due to background isotropic 
radiation pressure. 
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